AD.  A 0 3 1 2 6 8 


SSS-R-76-2865 


TELESEISMIC  COUPLING  FROM  THE  SIMULTANEOUS 
DETONATION  OF  AN  ARRAY  OF  NUCLEAR  EXPLOSIONS 

J.  T.  Cherry 
T.  C.  Bache 
W.  0.  Wray 
J.  F.  Masso 


Topical  Report 


Sponsored  by 

Advanced  Research  Projects  Agency 
ARPA  Order  No.  2551 


This  research  was  supported  by  the  Advanced 
Research  Projects  Agency  of  the  Department  of 
Defense  and  was  monitored  by  AFTAC/VSC,  Patrick 
AFB  FL  32925,  under  Contract  No.  F08606-75-C- 
0045. 

The  views  and  conclusions  contained  in  this  docu- 
ment are  those  of  the  authors  and  should  not  be 
interpreted  as  necessarily  representing  the  of- 
ficial policies,  either  expressed  or  implied,  of 
the  Advanced  Research  Projects  Agency,  the  Air 
Force  Technical  Applications  Center,  or  the  U.S. 
Government. 


Approved  for  Public  Release;  Distribution  Unlimited 


February  10,  1976 

P O BOX  1620,  LA  JOLLA.  CALIFORNIA  92038.  TELEPHONE  (714)  453  0060 


SYSTEMS,  SCIENCE  AND  SOFTWARE 


SSS-R-76-2865 


TELESEISM I C COUPLING  FROM  THE  SIMULTANEOUS 
DETONATION  OF  AN  ARRAY  OF  NUCLEAR  EXPLOSIONS 

J.  T.  Cherry 
T.  C.  Bache 
W.  0.  Wray 
J.  F.  Masso 


Topical  Report 


Sponsored  by 

Advanced  Research  Projects  Agency 
ARPA  Order  No.  2551 


This  research  was  supported  by  the  Advanced 
Research  Projects  Agency  of  the  Department  of 
Defense  and  was  monitored  by  AFTAC/VSC,  Patrick 
AFB  FL  32925,  under  Contract  No.  F08606-75-C- 
0045. 

The  views  and  conclusions  contained  in  this  docu- 
ment are  those  of  the  authors  and  should  not  be 
interpreted  as  necessarily  representing  the  of- 
ficial policies,  either  expressed  or  implied,  of 
the  Advanced  Research  Projects  Agency,  the  Air 
Force  Technical  Applications  Center,  or  the  U.S. 
Government. 


Approved  for  Public  Release;  Distribution  Unlimited 

February  10,  1976 


P.  O.  BOX  1620,  LA  JOLLA,  CALIFORNIA  92038,  TELEPHONE  (714)  453-0060 


UNCLASSIFIED 


5ECU  HI  TV  CLASS'FICATION  OF  This  "AGE  /i««  dnt*r,d) 


HEAD  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


REPORT  DOCUMENTATION  PAGE 


I.  REPORT  NIjMOER  l.  GOVT  ACCESSION  NO.  3.  PECl°TNT’S  CATALOG  NUMBER 


4.  TITLE  (»>d  Submit) 


TZLESEISMIC  .COUPLING  FRCM  THE  SIMULTANEOUS  \ 
/DETCNATICN  OF  AN  ARRAY  OF  NUCLEAR  EXPLOSIONS, 


8.  CONTRACT  OR  GRANT  NUMBERn) 


7.  vAuTMO^r*)  * 

J.  Ty/cherry,  T.  C./4ache,  W.  0./toray4m£-  '7 
J.  F^/Masso  IkJ_z=> 


f.  PERFORMING  organization  name  ano  aooress 

Systems,  Science  and  Software 

P.  0.  Bax  1620  ( 

La  Jolla,  California  92038 


II.  CONTROLLING  OFFICE  NAMe  ANO  AOORESS  S 7 

VELA  SEISMOLOGICAL  CENTER  ( (l 

312  Montgomery  Street  • 

Alexandria,  Virginia  22314 


IS.  MONITORING  AGENCY  name  » AOO  RESS<II  dllltrmnt  from  Controlling  OHIct)  IS.  SECURITY  CLASS,  (at  rAl*  rmport) 


13.  NUMBER  O'  RAGES 
26 


( \ y i 


Unclassified 


I5«.  OECL  ASSIFIC  ATI  ON/ OOWN  GRADING 
SCHEDULE 


18.  DISTRIBUTION  STATEMENT  (of  thie  Report) 


Approved  for  public  release;  distribution  unlimited 


17.  DISTRIBUTION  STATEMENT  (of  the  • betrect  entered  In  Block  20,  it  different  from  Report) 


19.  KEY  WOROS  (Continue  on  teeeree  ride  if  neceeeery  end  identify  by  block  number) 


Explosion  Seismology 
Multiple  Explosions 
Constitutive  Modeling 
Shock  Interaction 


Teleseismic  Magnitudes 
Earth  Structure 
Seismic  Coupling 


\ 


DO  | jan *73  1473  EOITION  OF  I NOV  «s  IS  OBSOLETE*" 

7-  • , ,c 


ified 


rity  CLASSIFICATION  of  THIS  r AOCXIFhon  Dot.  Ent. r.o 


-UNCLASSIFIED 

3gCj*tT>  CLASSIFICATION  OF  THIS  P AOCrWin  P>K  Enffd) 

r>  Enhanced  tension  failure  did  occur  in  the  multiple  shot  calculation. 
However,  for  this  to  seriously  affect  teleseismic  ground  motion  it  would 
require  an  explosive  spacing,  yield  and  material  properties  combination  in 
which  fracturing  from  the  array  is  rmch  greater  than  that  produced  by  a 
single  charge  of  equivalent  yield. -y 

'An  array  in  which  the  important  parameters  are  optimized  to  significantly 
alter  the  teleseismic  ground  motion  could  probably  be  constructed.  However, 
it  is  much  more  likely  that  anomalous  observations  of  teleseisnic  ground 
motion  from  an  explosive  array  cure  due  to  overburden  effects  and  high  (low) 
coupling  material  properties.  If  these  factors  fail  to  explain  the  anomaly, 
shock  interaction  calculations  should  be  performed  but  should  be  based  on  a 
good  estimate  of  material  strength  for  the  rock  at  the  test  site. 


l ACCS  WOK  lor 

•in 

mi"  Sertlu.i 

ODC 

loll  StClMI 

□ 

tmo»N«mcrt 

□ 

JiisiifiuriM  ..... 

IT  

(ISTIIIITIN/miUllliTT  CODES 

am.  mu. 

j 

a 

Ul 

& 

*1 

(V 


ii 


UNCLASSIFIED 


sicunity  classification  of  this  faok»ti«»  d«« 


R-2865 


TABLE  OF  CONTENTS 

Page 

I.  INTRODUCTION  1 

II.  NEAR  SOURCE  GROUND  MOTION  2 

III.  TELESEISMIC  GROUND  MOTION  11 

IV.  SUMMARY  AND  CONCLUSIONS 22 

APPENDIX  A 23 

REFERENCES 26 


ill 


R-2865 


I.  INTRODUCTION 

This  report  presents  the  results  of  a theoretical 
study  undertaken  to  determine  if  teleseismic  ground  motion 
may  be  enhanced  by  the  simultaneous  detonation  of  a closely 
spaced  array  of  nuclear  explosions.  The  specific  array 
analyzed  consisted  of  three  15  kt  explosions  equally  spaced 
165  meters  apart.  The  explosive  array  was  assumed  contained, 
i.e.,  coupling  effects  due  to  cratering  were  not  included 
in  the  analysis. 

For  this  specific  explosive  array  and  the  material 
properties  assumed  for  the  near  source  rock  environment  we 
found  no  significant  change  in  the  teleseismic  signature  of 
the  multiple  shot  array  over  that  produced  by  a single  explo- 
sive of  equivalent  yield.  This  was  true  in  spite  of  the  fact 
that  the  near  source  fracture  pattern  from  the  array  differed 
substantially  from  that  produced  by  a single  explosion. 

A general  answer  to  the  question  of  the  effect  of 
shock  interaction  on  teleseismic  ground  motion  is  not  avail- 
able from  our  study.  However,  our  results  indicate  that  if 
anomalous  teleseismic  magnitudes  are  observed  from  a particu- 
lar explosive  array,  then  the  explanation  should  first  be 
sought  in  the  material  properties  and  the  overburden  pressure 
associated  with  the  specific  rock  type  and  depth  of  burial. 

It  may  be  possible  for  shock  interaction  to  contribute 
to  anomalous  teleseismic  magnitudes  for  an  explosive  array. 
However,  for  this  to  occur  tension  fracturing  for  the  array 
must  be  much  greater  than  that  from  a single  explosion 
of  equivalent  yield.  Since  material  strength  is  the  key 
parameter  which  determines  the  extent  of  fracturing,  it  is 
very  important  that  a firm  estimate  of  this  parameter  be  avail- 
able before  attempting  to  evaluate  the  effect  of  the  shock 
interaction. 
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II.  NEAR  SOURCE  GROUND  MOTION 


Details  of  the  techniques  used  to  predict  the  tele- 
seismic  ground  motion  from  an  explosive  array  have  been  pre- 
sented previously  (Cherry,  et  ah,  1975a;  Bache,  et  al. , 1975). 
Basically  the  technique  requires  that  the  divergence,  V • S, 
and  curl,  V x S,  of  the  displacement  field  be  monitored  over 
a spherical  surface  which  is  outside  the  nonlinear  material 
response  region.  The  time  histories  of  V • S and  V x S are 
then  decomposed  into  a multipole  expansion  in  spherical 
harmonics.  This  expansion  provides  an  equivalent  point 
source  representation  which  may  be  used  to  quantitatively 
establish  the  seismic  coupling  of  the  array. 

Figure  2.1  shows  a sketch  of  the  explosive  array, 
the  nonlinear  region  and  the  spherical  surface  over  which 
V • S and  V x s were  monitored.  The  computational  grid  ex- 
tended far  enough  beyond  this  surface  so  that  external 
reflections  would  not  obviate  the  free  field  assumption. 

The  tensile  fracture  pattern  around  two  adjacent  cavi- 
ties is  shown  in  Figures  2.2  through  2.5.  Note  the  "pre- 
splitting" between  the  cavities  and  the  alteration  of  the 
fracture  pattern  from  that  expected  from  a spherically  sym- 
metric explosive  source.  This  last  item  is  illustrated 
dramatically  in  Figure  2.6  which  shows  the  distribution  of 
fractures  in  the  entire  grid  at  11.57  msec.  The  portion  of 
the  grid  displayed  in  Figure  2.5  is  a subset  of  that  shown 
in  Figure  2.6.  The  fracture  pattern  above  and  to  the  left 
of  the  outer  source  is  that  expected  from  a single  explosive. 
The  alteration  of  this  fracture  pattern  is  caused  by 
interaction  of  the  stress  fields  from  adjacent  explosives. 

The  extent  of  fracturing  is  controlled  by  the  assumed 
material  strength.  For  this  calculation  we  assume  that 
material  strength  (Y)  was  related  to  stress  state  (P)  by 


2 


Spherical  surface  over 


Sketch  of  the  explosive  array,  the  resulting 
nonlinear  region  and  the  surface  over  which 
the  properties  of  the  outgoing  displacement 
field  were  monitored. 
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Figure  2.2. 


Crack  location  and  orientation  4.93  msec  after 
detonation.  The  lines  indicate  the  crack 
orientation  in  the  plane  while  a dot  indicates 
a zone  in  which  tensile  failure  out  of  the  plane 
has  occurred. 
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Figure  2.3.  Crack  location  and  orientation  6.85  msec  after 
detonation. 
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Figure  2.6.  Complete  crack  location  and  orientation  11.57 
msec  after  detonation. 
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" Y,  + Ym  ' P > Pm  ' U-U 

= 0,  Y < 0 , 


where 


Y = 0.4  kbar, 

o ' 

Ym  = 0.15  kbar, 

P =0.45  kbar . 

Ill 

The  justification  for  assuming  Y = f (P)  along  with 
the  tension  failure  model  used  when  material  strength  is 
exceeded  has  been  given  by  Cherry,  et  al.  (1975b) . For  our 
purposes  it  will  be  sufficient  to  point  out  that  a good  mea- 
sure of  the  tensile  strength  is  given  by  Y /2;  i.e.,  a ten- 

0 

sile  strength  of  0.2  kbar  was  assumed  for  this  calculation. 

Our  main  objective  was  to  determine  if  the  fracture 
pattern  unique  to  a multiple  shot  array  could  significantly 
effect  teleseismic  ground  motion.  The  results  shown  in 
Figure  2.6  suggest  that  in  the  near  field  the  multiple  shot 
is  quite  different  from  a single  explosion.  In  the  next 
section  it  will  be  seen  how  these  differences  are  transmitted 
to  the  far  field.  However,  we  must  point  out  that  the  sensi- 
tivity of  our  results  to  tensile  strength  remains  to  be  demon- 
strated. 

Finally,  we  note  that  the  maximum  allowable  stress 
difference  assumed  in  this  calculation  was 

Y + Y =0.55  kbar. 
o m 
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This  value  is  identical  to  that  used  for  a previous  calcula- 
tion, reported  by  Cherry,  et  al . (1975c).  However,  no  tension 

failure  was  allowed  in  the  previous  calculation.  Taking  the 
calculation  described  here  together  with  the  one  reported  by 
Cherry,  et  al . (1975c)  , we  have  two  sources  which  will  be 
described  as  the  "no  tensile  failure"  source  and  the  "with 
tensile  failure"  source.  In  the  first  the  multiple  explosion 
array  is  different  from  a single  explosion  only  in  the  asym- 
metry of  the  source  arrangement.  In  the  calculation  described 
here,  the  effect  of  tensile  failure  is  included. 

In  our  previous  report  (Cherry,  et  al. , 1975c)  we 
showed  that  the  "no  tensile  failure"  multiple  source  calcula- 
tion led  to  a teleseismic  signature  that  was  little  different 
from  that  for  a single  explosion.  In  the  following  section 
the  seismic  radiation  for  both  multiple  sources  will  be  dis- 
cussed and  compared  to  determine  if  teleseismic  ground  motion 
is  significantly  altered  by  the  multiple  shot  array. 
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III.  TELESEISMIC  GROUND  MOTION 


In  order  to  study  the  far  field  seismic  signature  of 
an  explosion  event  like  that  considered  here,  it  is  necessary 
to  merge  the  nonlinear  finite  difference  source  calculations 
presented  in  Section  II  with  the  elastic  wave  propagation 
methods  of  theoretical  seismology.  This  is  done  via  an 
equivalent  elastic  source  representation  of  the  multiple 
explosion  source. 

The  familiar  reduced  displacement  potential,  RDP , 
representation  of  a spherically  symmetric  explosion  is  an 
elementary  equivalent  elastic  source.  For  more  complex 
sources  a representation  in  terms  of  an  expansion  of  the  out- 
going displacement  field  in  spherical  harmonics  is  appropriate. 
The  calculation  of  an  equivalent  elastic  source  in  this  latter 
form  for  a complex  asymmetric  explosion  source  was  first  ac- 
complished by  Cherry,  et  al . (1975a) . A detailed  presentation 
of  the  theory  and  its  implementation  is  given  by  Bache,  et 
al . (1975).  For  convenience,  a summary  is  given  in  Appendix 

A. 

Using  Eq.  (A. 5)  of  the  appendix,  an  equivalent  elastic 

source  was  computed  for  the  multiple  explosion  source.  The 

. (4 ) 

nonzero  multipole  coefficients  included  the  monopole,  A (u) , 

o o 

which  is  the  spherically  symmetric  portion  of  the  source. 

The  nonzero  higher  order  terms  were 

} (w),  (u)  = -A^  (o>),  £ = 2,  4,  6,  ...  , 

Where  the  3-axis  was  taken  to  be  along  the  line  of  explosion 
centers.  The  terms  with  £ = 2 are  referred  to  as  the  quad- 
rupole  terms  and  those  with  £ = 4 as  the  octupole  terms, 
etc. 

The  main  objective  is  to  determine  the  teleseismic 
signature  of  the  multiple  explosion.  Therefore,  it  is 
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appropriate  to  draw  conclusions  from  the  far  field  portion 
of  the  displacement  field;  that  is,  only  terms  decaying  with 
R-1  need  be  retained  in  the  expansion,  Eq.  (A. 5).  In  this 
case  the  P and  S wave  displacements  reduce  to 


u (R,oj) 


-ik  R 
e P_ 

k2R 

P 


A(4)  (u>) 
2 0 


(3cos  28+1)  + H_0_T 


u (R,oj) 
s ■* 


2e 


-ik  R 
s 


k2R 


I-  1 

\ 2 


B (1)  (ui) 
2 1 


sin  20  + H.O.T, 


S' 


(3.1) 


where  H.O.T.  indicates  the  higher  order  terms  in  the  infinite 
series.  In  the  frequency  range  of  interest,  terms  higher 
than  the  quadrupole  have  almost  no  influence.  The  geometry 
is  shown  in  Figure  3.1.  With  this  geometry  the  displace- 
ments are  independent  of  the  "azimuthal"  coordinate,  <t> . 

The  most  direct  way  to  study  the  equivalent  elastic 
source  is  to  examine  the  multipole  coefficients  themselves. 

The  amplitude  spectral  density  of  the  monopole  and  quadru- 
pole coefficients  is  plotted  in  Figure  3.2.  The  coefficients 
for  both  source-calculations,  with  and  without  tension  failure, 
are  shown.  From  this  plot  the  following  is  apparent: 

• The  monopole  differences  are  mainly  at  short 
periods  where  the  tensile  failure  model  is 
larger  (e.g.,  by  a factor  of  1.7  at  1 Hz). 

(4 ) 

• The  P wave  quadrupole  (Aj0  ) is  an  order  of 
magnitude  less  than  the  monopole  for  fre- 
quencies less  than  1 Hz  for  both  sources. 

For  higher  frequencies  the  P wave  quadrupole 
is  increasingly  more  important.  Selected 
ratios  of  |a|$ ) | to  |Aj*) | are  as  follows: 


With  Tension 
Failure 

No  Tension 
Failure 

2 Hz 

1.8 

4.3 

3 Hz 

0.7 

2.0 

4 Hz 

0.6 

1.2 

12 

Amplitude  Spectral  Density 


e 3.2.  The  amplitude  spectral  density  of  the  multipole 
coefficients  (monopole  and  quadrupole)  for  two 
multiple  burst  calculations,  with  and  without 
tensile  failure. 
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• The  S wave  quadrupole  (B2 1 ) is  relatively  constant 
with  respect  to  the  monopole.  Selected  ratios  of 
| A ^4 ) | to  |b  (1)  | are  tabulated  below: 

With  Tension  No  Tension 


Failure  Failure 

0.05  Hz  3.8  6.8 

0.5  Hz  3.1  6.2 

1.0  Hz  2.4  5.8 

2.0  Hz  2.7  6.3 


From  the  tabulated  values  presented  it  is  clear  that 
tension  failure  does  enhance  the  quadrupole  perturbation  on 
the  spherically  symmetric  portion  of  the  source.  The  enhance- 
ment is  by  a factor  of  two  or  so.  However,  except  at  high 
frequencies,  the  quadrupole  terms  remain  too  small  to 
have  much  effect  on  the  teleseismic  signal. 

A number  of  radiation  patterns  for  the  two  sources 
are  shown  in  Figures  3.3  and  3.4.  The  displacement  spectral 
components  in  these  patterns  are  normalized  to  those  due 
to  the  monopole  alone.  The  higher  order  perturbation  on  the 
spherically  symmetric  source  is  thereby  directly  displayed 
and  is  independent  of  R. 

Examining  the  radiation  patterns  in  Figure  3.3,  it 
is  clear  that  the  "no  tensile  failure"  source  behaves,  in 
terms  of  the  far  field  radiation,  almost  the  same  as  a 
spherically  symmetric  source.  As  was  seen  from  the  multi- 
pole coefficients,  the  addition  of  tensile  failure  increases 
the  quadrupole  component  but  it  is  still  quite  small  (Figure 
3.4).  Even  at  2 Hz  the  radiated  P wave  for  the  "with  ten- 
sile failure"  source  is  no  more  than  15  percent  different 
from  that  due  to  the  monopole  alone.  The  S waves  are,  of 
course,  much  too  small  to  contribute  a significant  sP  phase 
to  the  early  arriving  body  wave  recording  from  which  m^  mea- 
surements are  made.  However,  it  is  true  that  rather  large 
(for  an  explosion)  direct  S waves  would  be  transmitted  to 
the  teleseismic  field. 
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(a)  0.5  Hz 


.1  '<•  T-« 


(b)  1.0  Hz 


Figure  3.3.  Radiation  patterns  for  the  "no  tensile  failure" 
multiple  source  calculation.  These  patterns  show  the  far 
field  P and  S wave  displacement  spectral  components  normalized 
to  that  from  the  monopole  contribution  above.  While  included 
in  the  calculation, terms  higher  than  the  quadrupole  make  al- 
most no  contribution  at  these  frequencies.  These  patterns 
are  symmetric  about  the  line  of  explosions  which  is  indicated 
as  the  z-axis. 
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From  these  radiation  patterns  and,  as  a double  check, 
from  a number  of  synthetic  seismogram  calculations,  it  was 
concluded  that  the  multiple  source  appears,  as  far  as  tele- 
seismic  P wave  and  Rayleigh  wave  recordings  are  concerned, 
like  a single  spherically  symmetric  explosion.  The  only  re- 
maining question  regards  the  yield  of  this  equivalent  single 
source.  This  can  be  determined  by  computing  the  RDP  from  the 
monopole  using  Eq.  (A. 8)  of  Appendix  A.  The  "equivalent 
RDP"  for  the  multiple  explosion  can  then  be  compared  to  that 
obtained  by  a linear  superposition  of  three  single  sources. 

RDP  calculations  for  the  tension  failure  and  no  tension  failure 
cases  are  shown  in  Figure  3.5.  The  tension  failure  RDP  shows 
the  "overshoot"  reported  by  Cherry,  et  al.  (1975b). 

In  Figure  3.6  and  3.7  the  amplitude  of  the  transformed 

A 

reduced  velocity  potential,  |4'(cj)|,  from  the  1-D  calculations 
(Figure  3.5)  is  compared  to  the  equivalent  quantity  for  the 
two  multiple  explosions  under  study.  The  linear  superposition 
is  done  by  (1)  scaling  the  source  as  if  it  were  a single  45 
kt  explosion  and  (2)  considering  the  event  to  be  the  super- 
position of  three  15  kt  explosions.  Due  to  the  cube  root 
scaling  laws,  the  only  difference  in  the  two  is  that  the 
second  representation  is  somewhat  richer  in  high  frequencies. 

From  Figures  3.6  and  3.7,  we  see  that  for  neither 
case  does  the  "equivalent  RDP"  from  the  spherically  symmetric 
portion  of  the  multiple  explosion  calculation  differ  sub- 
stantially from  that  obtained  by  a linear  superposition  of 
one-dimensional  sources.  It  is  interesting  to  note  that  the 
corner  frequency  for  the  no  tensile  failure  calculation  falls 
between  that  for  the  3 * 15  kt  and  45  kt  superposition 
sources.  For  the  second  calculation,  with  tensile  failure, 
the  "equivalent  RDP"  is  very  little  different  than  that  for 
the  45  kt  source. 
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Figure  3.5.  The  amplitude  of  the  transformed  reduced  velocity 
potential  for  two  spherically  symmetric  explosion 
calculations.  The  explosion  yield  is  15  kt. 
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Figure  3.6. 


For  the  no  tensile  failure  calculation,  the 
"equivalent  RDP"  representation  of  the  multiple 
shot  is  compared  to  the  source  obtained  by  a linear 
superposition  of  spherically  symmetric  explosions. 
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Figure  3.7.  Comparison  of  the  "equivalent  RDP" 
tions  of  the  multiple  explosion  foj 
culations  including  the  effect  of  1 
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IV.  SUMMARY  AND  CONCLUSIONS 

Calculations  performed  to  date  have  produced  no 
significant  enhancement  of  teleseismic  ground  motion  from  a 
multiple  shot  a* ray.  This  was  true  in  spite  of  the  fact 
that  the  near  source  fracture  pattern  from  the  array  dif- 
fered significantly  from  that  produced  by  a spherically 
symmetric  explosion  source.  Fairly  large  perturbations  in 
the  seismic  radiation  pattern  did  appear  but  at  frequencies 
outside  the  teleseismic  band. 

Enhanced  tension  failure  did  occur  in  the  multiple 
shot  calculation.  However,  for  this  to  seriously  affect 
teleseismic  ground  motion  it  would  require  an  explosive 
spacing,  yield  and  material  properties  combination  in  which 
fracturing  from  the  array  is  much  greater  than  that  produced 
by  a single  charge  of  equivalent  yield. 

An  array  in  which  the  important  parameters  are  opti- 
mized to  significantly  alter  the  teleseismic  ground  motion 
could  probably  be  constructed.  However,  it  is  much  more 
likely  that  anomalous  observations  of  teleseismic  ground 
motion  from  an  explosive  array  are  due  to  overburden  effects 
and  high  (low)  coupling  material  properties.  If  these  fac- 
tors fail  to  explain  the  anomaly,  shock  interaction  calcula- 
tions should  be  performed  but  should  be  based  on  a good  esti- 
mate of  material  strength  for  the  rock  at  the  test  site. 
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APPENDIX  A 

EQUIVALENT  ELASTIC  SOURCE  REPRESENTATION 


The  Fourier  transformed  equations  of  motion  in  a homo- 
geneous, isotropic,  linearly  elastic  medium  may  be  written 


where  u is  particle  displacement  and  k and  kQ  are  the 
~ a p 

compressional  and  shear  wave  numbers.  The  Cartesian  poten- 
tials x(4)  and  x are  defined  by 

X(4)  = V • u 

(A.2) 

X = \ V x u 


and  may  be  easily  shown  to  satisfy  the  wave  equation 

V2-(j)  + k2^-(j)  = o , j = 1,  2,  3,  4,  (A. 3) 

where  k = k = u/a  , k-  = k0  “ w/B  , i = 1,2,3.  This  equa 

4 U IP 

tion  has  as  a solution  the  following  expansion  in  spherical 
eigenfunctions  (e.g.,  Morse  and  Feshbach  (1953)), 


00 


E hP)(>'aR) 


cos  m$ 


(A.  4) 


♦ 


C“> 


sin  m<J> 


Pj(cos0) 


where  the  h^2^  are  spherical  Hankel  functions  of  the  second 
kind  and  the  P™  are  associated  Legendre  functions.  The 
vector  R has  as  components  the  spherical  coordinates  R,0,<{>. 


Equations  (A. 4) together  with  (A.l)  provide  an  equiva- 
lent elastic  source  representation  of  the  outgoing  elastic 
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waves.  The  values  of  the  multipole  coefficients  (u>)  , 

B^jj^  (oj)  , j = 1,2, 3,4,  prescribe  the  displacement  field  at  all 
points  in  the  homogeneous  elastic  medium  where  (A.l)  applies. 

The  point  source  representation,  Eq.  (A. 4)  , provides  a 
generalization  of  the  commonly  used  center  of  dilatation 
(i.  = 0),  couple  (£  = 1)  and  double  couple  (1  = 2)  point 
sources,  as  well  as  more  complex  sources.  For  example,  for 
a horizontal  double  couple  source  of  unit  amplitude,  the  non- 
zero coefficients  in  (A. 4)  are 

-A^  = B^  « A^  = (B2/ a2)  B^  . 

21  21  22  J 22 


Given  the  displacement,  u ^ , or,  alternatively,  the 
potentials  x ^ » one  may  determine  the  multipole  coefficients 
using  (A. 4)  • For  example,  using  the  orthogonality  of  the 
spherical  eigenfunctions,  one  may  derive 


*ii5  <“> 


CkaR> 


I J X^CR.w)  (cos  0 ) < >sin0d0d4> 


cos  m<J> 


sm  m<j> 


(A.S) 


where 


’ll*!) (H-m) ! 
2 it  (l+m) 


, m ^ 0 , 


(2i+l)/4ir  . 


That  is,  given  the  displacements,  u(R,w)  on  a sphere  of 
radius  R,  one  can  use  (A. 2),  (A. 5)  and  (A. 4)  to  compute  the 
displacements  at  any  point  in  the  elastic  regime.  This  pro- 
cedure was  first  suggested  for  geophysical  applications  by 
Archambeau  (1968). 

Using  the  equivalent  elastic  source  representation, 
highly  nonlinear  finite  difference  calculations  of  the  near 
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source  earthquake  or  explosion  displacement  fields  may  be 
coupled  to  the  elastic  field  where  efficient  elastic  wave 
propagation  techniques  can  be  applied.  Carrying  the  non- 
linear computation  out  to  a radius  at  which  the  material 
response  is  in  the  elastic  regime,  Eq.  (A. 5)  may  be  applied 
to  compute  the  multipole  coefficients  and,  thereby,  the 
equivalent  elastic  source. 

For  a spherically  symmetric  explosion,  the  explosion 
induced  elastic  field  is  often  represented  by  the  reduced 
displacement  potential,  ip(r),  defined  by 


ur(r,t) 


ICll  + 

R2 


1 

R~ 


dr 


(A. 6) 


where  t = t - R/a,  is  the  retarded  time.  Note  that  ¥(t) 
is  independent  of  R. 

Applying  the  Fourier  transform,  (A. 6)  becomes 


(1  + ik  r)  -ik  R _ 

u_(R,uO  « ~ — e a 

A 2 


(A. 7) 


Then,  carrying  out  the  algebra  (observing  that  the  u of  (A. 2) 
is  the  Cartesian  displacement  vector),  it  may  be  shown  that 
the  multipole  coefficients  reduce  to  a single  nonzero  term, 
the  monopole: 

A^C)  ■ -ik*?((D)  . (A. 8) 

0 0 H 
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